NUSc/93-04 

Van Hove Exciton-Cageons and High-Tc Superconductivity: XA: 
Role of Spin-Orbit Coupling in Generating a Diabolical Point 

R.S. Markiewicz 

Physics Department and Barnett Institute, Northeastern U., Boston MA 02115 

Abstract 

Spin-orbit coupling plays a large role in stabilizing the low-temperature or- 

thorhombic phase of La2-xSra;Cu04. It splits the degeneracy of the van Hove 

singularities (thereby stabilizing the distorted phase) and completely changes 

the shape of the Fermi surfaces, potentially introducing diabolical points into 

the band structure. The present paper gives a detailed account of the resulting 

electronic structure. 

A slave boson calculation shows how these results are modified in the 

presence of strong correlation effects. A scaling regime, found very close to the 

metal-insulator transition, allows an analytical determination of the crossover, 

in the limit of zero oxygen-oxygen hopping, too 0. Extreme care must 

exercised in chosing the parameters of the three-band model. In particular, 

too is renormalized from its LDA value. Furthermore, it is suggested that 

the slave boson model be spin-corrected, in which case the system is close to 

a metal-insulator transition at half filling. 
PACS numbers : 71.27.+a, 74.20.Mn 
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I. INTRODUCTION 



The suggestion that the van Hove singularity (vHs) plays an important role in high-T^ 
superconductivity has been greatly strengthened by recent photoemission studies, which 
find that the Fermi level is very close to the vHs at the compositions of optimum Tc, both 
in BiaSrsCaCusOg (Bi-2212) § and in YBasCugOr (YBCO) §. Furthermore, it has been 
shown theoretically that the vHs's in La2_xSr2:Cu04 (LSCO) and related compounds drive 
a very strong electron-phonon interaction The mechanism is via a vHs-Jahn- Teller 

effect, wherein the two vHs's provide the electronic degeneracy which can be split by a 
lattice distortion - a band Jahn- Teller (JT) effect. However, in order to reproduce the 
experimental phase diagram, it was found that electronic correlations play an essential 
role, by renormalizing the band structure. The present paper is a first attempt to treat both 
the electron-hole coupling and correlation effects on an equal footing. 

It is found that spin-orbit coupling plays a dominant role in driving the structural insta- 
bility, and leads to major reconstruction of the Fermi surfaces. It acts like a form of umklapp 
scattering, splitting the degeneracy of the vHs's (even in the LTO phase!), and producing 
gaps in the vicinity of the old vHs's, thereby transforming the large Fermi surfaces into small 
pockets in better agreement with Hall effect measurements. However, while this coupling 
can restructure the Fermi surfaces, it cannot eliminate vHs's. Instead, the vHs's are split 
and shifted in energy, forming a very complex pattern of up to four vHs's, with correspond- 
ing peaks in the density of states (dos). A detailed description of the various possibilities is 
presented. 

A most interesting possibility is that, if spin-orbit coupling is strong enough, the band 
may posess a diabolical point, at which the Fermi surface shrinks to a point, and the dos 
vanishes. The possible role of such a state in the antiferromagnetic insulating phase at half 
filling is discussed. 

The importance of these features depends critically on relative values of several band 
parameters - both bare parameter values and how the values renormalize in the presence of 
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correlation effects. It is found that there remains considerable uncertainty in the optimal 
values for several parameters, but approximate estimates are given. 

Finally, correlation effects are included by means of a slave boson calculation. In the 
absence of spin-orbit (and other electron-phonon) coupling, a full calculation is given (to 
lowest order in 1/N) for the three-band model of the Cu02 planes, and the metal-insulator 
phase boundary is found. In the limit of strong renormalization of the Cu-0 hopping pa- 
rameter, a scaling regime is found in which the self-consistent equations greatly simplify. 
The phase boundary is found to lie in the range of the optimum parameters. 

When spin-orbit and electron-phonon coupling are included, it is found that there is 
still a scaling regime, and that, since the additional coupling is small, the phase boundary 
is nearly unchanged. For the parameters chosen, the state at half filling is semimetallic; 
however, a larger, renormalized value of the spin-orbit coupling parameter might be more 
appropriate. In this case, the state at half filling would become a zero-gap semiconductor 
with diabolical point. 

This paper naturally divides into two parts. First, Sections II, III show how spin-orbit 
(and electron-phonon) coupling modifies the three-band, tight-binding model of the Cu02 
planes, splitting the degeneracy of the vHs's even in the LTO phase. The new Fermi surfaces 
and vHs's are described in detail, and it is shown that at half filling there can be either 
compensated, semimetallic Fermi surfaces, or a zero-gap semiconductor with a diabolical 
point at precisely half filling (Section III). The remainder of the paper introduces correlations 
into the problem, via slave boson calculations. This includes a discussion of the proper 
choice of parameters (Section IV), the self-consistent equations and their modification in 
the presence of spin-orbit coupling (Section VA), the analysis of a scaling regime (and its 
modification in the presence of spin-orbit coupling) which exists near the metal-insulator 
transition (Section VB,C), and a more general discussion of parameter renormalization near 
this transition (Section VD). Conclusions are discussed in Section VI, while Section VII 
discusses a number of new features which arise in the presence of spin-orbit coupling and 
the possible diabolical point. 
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II. ROLE OF SPIN-ORBIT COUPLING IN THE STRUCTURAL TRANSITIONS 



Spin-orbit coupling plays two major roles in the structural transitions. First, since 
spin-orbit coupling is linear in the tilt angle 9, while ordinary electron-phonon coupling 
is quadratic, the spin-orbit term dominates in driving the structural transition Sec- 
ondly, spin-orbit coupling can split the degeneracy of the vHs's, leading to a large electronic 
contribution to the stabilization of the LTO phase [^]. 



A. Spin-Orbit Coupling 

In La2_a;Ba3;Cu04 (LBCO), there is a series of structural transitions as the temperature 
is lowered, from a high-temperature tetragonal (HTT) phase near room temperature, to the 
low-temperature orthorhombic (LTO) phase, and, in a restricted doping range, to the low- 
temperature tetragonal phase (LTT), which has a deleterious effect on superconductivity. It 
has long been recognized that the distortion in the LTT phase splits the degeneracy of the 
two vHs's [H, thereby acting like a classical charge density wave, or static JT effect. 



On the other hand, it was generally assumed [10| that since the two vHs's remain degen- 
erate by symmetry in the LTO phase, the LTO transition is unrelated to the vHs's. This 
argument is incorrect, and in fact there are two distinct ways in which splitting of the vHs's 
can drive the LTO transition [Q. From general symmetry arguments, it can be shown that 
spin-orbit coupling can lift the orbital degeneracy on most of the orthorhombic Brillouin 
zone, leading to a splitting of the vHs density of states peak. Alternatively, the LTO phase 
could be a dynamic JT phase, in which the local symmetry differs from the global one, 
stabilized by a dynamic splitting of the vHs degeneracy. Based on this picture, I presented 
0] a simplified calculation of the resulting doping dependence of the LTO phase transition 
and the competing superconductivity. 

A complication in analyzing this transition is that the soft mode of the transition in- 
volves a tilting of the CuOg octahedra, which has conventionally been found to couple only 
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quadratically to the electrons (non-Migdal behavior) [p],p|,p!2|. Recently, however, it has 
been pointed out that spin-orbit coupling also introduces a linear coupling between the elec- 
trons and the tilt mode [Q. In the present paper, I explore the role of spin-orbit coupling in 
driving the LTO transition. This Subection describes the modifications of the band struc- 
ture due to spin-orbit coupling in the LTO and LTT phases. Unlike earlier work no 
magnetic ordering is assumed. 

The analysis is based on the standard three-band model of the Cu02 plane, involving 
both Cu-0 (tcuo) and 0-0 (too) hopping: 

H = m^Adj + E^c«o[4p,+5 + (c-c.)] + E^oo[p]+^P,+5. + (c.c.)] 

+Unj^nji), (1) 

where d) (p^) is a creation operator for holes on Cu (O), j is summed over lattice sites, 5 
over nearest neighbors, 5' over next-nearest (0-0) neighbors, and c.c. stands for complex 
conjugate. Energies are measured from the center of the O bands, with A the bare Cu-0 
splitting, and U is the on-site Coulomb repulsion. In a slave boson calculation, U produces 
correlations which renormalize the one-electron band parameters. Hence, in the present 
section the Coulomb term will be neglected, in order to study how the band structure is 
modified as extra terms are added to Eq. 1. In Section IV, the slave boson technique will 
be used to determine the renormalized band structure. 

Coffey, et al. showed that, in the presence of octahedral tilt, spin-orbit coupling 
introduces an additional, spin-dependent hopping term into the Hamiltonian, Eq. 1, of the 
form 

H' = Y. T.idU\j+5 ■ <PP,+lp + (c-c.)), (2) 

where a and (3 are spin indices, (Tq,^/3 is the vector of spin matrices, and the coupling parameter 
has the form 

4,+£ = (-l)^'^+^'-(Ai,A^,0), (3a) 
A,.+s = -(-1)^'^+^'^(A2,A;,0), (3&) 
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where the Aj's are microscopic parameters, estimated to be Xi = 'jiO, with 6 the octahedral 
tih angle and 72 > 7i > 0, = Aj in the LTO phase, and = in the LTT phase, and the 
values of Aj are larger by a factor in the LTT phase, for a given tilt angle. 

B. Modified Band Structure in the Presence of Spin-Orbit Scattering 

The LTO phase of LSCO has symmetry group Bmab; in the absence of spin-orbit cou- 
pling, this means that the electronic states on the X* face of the Brillouin zone have an 
additional two-fold degeneracy (beyond the usual spin degeneracy), which means that the 
degeneracy of the two vHs's cannot be split by a transition to a uniform orthorhombic phase. 

Spin-orbit scattering dramatically alters this situation. The spin ceases to be a good 
quantum number, and the eigenfunctions mix up and down spin states. There is still a 
two-fold Kramers degeneracy, but the extra degeneracy on the X* face is lifted everywhere, 
except at X* itself [0. This section will provide a detailed calculation of how the Fermi 
surfaces are changed by the spin-orbit coupling, Eq. 2. 

The orthorhombic unit cell is double the area of the HTT cell, containing 2 Cu's and 4 
O's per layer. Since up and down spins must be treated separately, the resulting Hamiltonian 
matrix is a 12 x 12. By properly symmetrizing the two Cu's (and corresponding O's) this 
can be reduced to two 6x6 matrices, which are Kramers' doublets. The eigenvalues are 
then given by 
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with t = tcuo, Si = sin{kia/2), Ci = cos{kia/2), i = x,y,u = AtooSxSy, v = itooCxCy, A± = 
A, and, for LTO, A^ = (Ai + iAa), A^, = -(A2 + iXi) (for LTT, A^ = y^Xi, Xy = -V2A2). 
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For the general case, Eq. 4 must be solved numerically, but in the special case too = 0, the 
determinant can be reduced to a quadratic equation 

[E{A -E)+ At^ + 2X^f = 4[t2(c^ + Cy) + 2\\^\'^c^ + 2\XyW]^ 

-IQt'c^CyX^ (5) 

with = + A|, A = Ai + A2, and Cj = cos{kia), i = x,y. When = Cy = 0, the 
right-hand side of Eq. 5 vanishes, leading to the (symmetry-required) two-fold degeneracy 
at X* and Y* (M, in the pseudo-tetragonal cell we are using). The generalization of Eq. 5 
to arbitrary too is given in Appendix I. 



C. Electron-Phonon Interaction 

As discussed elsewhere |Q, ordinary electron-phonon coupling can split the vHs degener- 
acy in the LTT phase, but not in a uniform LTO phase. Moreover, the coupling is quadratic 
in the tilt angle, 6, and hence of non-Migdal form However, after the spin-orbit 

coupling leads to a finite tilt, this ordinary electron-phonon coupling can enhance the tilt 
instability. Hence, in this subsection, the above expressions are modified to include a JT-like 
phonon coupling. In Ref. 0, I showed that an important contribution comes from coupling 
between the Cu dx2-y2 orbital and the O pz orbital. This contribution can be incorporated 
into the interaction matrix, Eq. 4, by the substitution (in the LTO phase) 

A± = A + r][2±{a, + Cy)], (6a) 

where 

ri = ^, (66) 

and t' (oc sin6) is the Cu-Op^ hopping matrix element. Note that strictly speaking, rj is 
inversely proportional to the energy, and hence will modify the eigenvalues of Eq. 4 at each 
angle. However, it should be sufficiently accurate for present purposes to treat as a fixed 



7 



parameter, of approximate magnitude t] ~ 2t'^/A. In this approximation, Equation 6 is 
very similar to the electron-phonon couphng recently proposed by Song and Annett |12[. 

From the form of Eq. 6, the electron-phonon coupling does not split the vHs degeneracy 
in the LTO phase - in fact, Eq. 6 reduces to a constant along the Brillouin zone diagonal, 
Cx + Cy = 0. In the LTT phase, the expression should be modified, as 

2 — (cj: + Cy) 2{sin^9xsl + sin'^OySy), 

where 6i is the tilt in the LTT phase (e.g., in one domain, 6y would vanish). 

D. Numerical Results 

Equation 4 has been solved numerically, and the resulting energy bands and Fermi sur- 
faces are illustrated in Figures 1-6. Figure 1 shows the resulting energy dispersions for a 
series of different band parameters. The appropriate choice of parameters for calculating the 
effects of correlations will be discussed in detail in Section IV. Here, for illustrative purposes, 
values close to the bare band parameters will be assumed. These can be taken as t = 1.3eV, 
A = AeV, too = 0.65eV^, A2 < Ai ^ 6meV, and 77 ~ 5 — lOmeV. The special points of 
the Brillouin zone are based on the tetragonal supercell, which has twice the area of the 
orthorhombic zone: F = (0, 0), X = (vr/a, 0), F = (0, 7r/a), M = (vr/a, vr/a), and M = M/2. 

The original, untilted case, with Aj = r/ = 0, is illustrated in Figure la. It can easily 
be understood as the usual dispersion found in the tetragonal cell (solid line) modified by 
the zone folding associated with the real space cell size doubling, which maps M ^ F. 
The corresponding Fermi surfaces are shown in Fig. 2a. The energy levels freely intersect 
without interacting along the orthorhombic zone boundary. In particular, they need not 
be perpendicular to this boundary when they intersect it. Under these circumstances, it is 
convenient to work in a double cell (the original tetragonal zone) and ignore the zone folded 
bands. This is what is usually done. 

However, spin-orbit coupling breaks this degeneracy almost everywhere along the zone 
boundary. This coupling acts like an umklapp term, mixing the original and zone folded 
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bands. The resulting, reconstructed Fermi surfaces bear little resemblance to the original 
surfaces in the regions of overlap. Figure 2b illustrates how the Fermi surfaces of Fig. 2a 
are modified by the presence of a small spin-orbit term (A2 = 6meV, Xi = rj = 0). The 
two lower energy levels {Ep = 4.5, 5eV) did not overlap, and are essentially unchanged. 
The third level, 5.28eV, was initially at the vHs, with area proportional to 1 + x (with x 
the Sr doping), but the overlap has distorted the Fermi surface into a pair of closed hole 
Fermi surfaces (recall that the figure illustrates only one-quarter of the tetragonal cell), each 
of area x/2. For the fourth energy level, 5.5eV, the reconstruction converts the original 
Fermi surface into an electron-like Fermi surface, centered at X and Y, and a hole-like 
surface, centered at M (and symmetry-equivalent points). The interaction has split the vHs 
degeneracy, as can be seen in Fig. lb, the energy dispersion corresponding to the Fermi 
surfaces of Fig. 2b. A gap has been opened up at the X-point of the zone, the location of 
the original vHs. However, the interactions cannot eliminate a vHs, but only split it into 
multiple vHs's. This will be further discussed in the following subsection. 

For convenience in discussing the effects of various parameters on the band structure, I 
have used large values of the parameters to enhance the effects (as will be discussed later, 
such larger values may indeed be more appropriate). Figure Ic shows that when A2 is 
increased to 50meV, the band structure is qualitatively unchanged, with only an enhanced 
band gap. As shown in Fig. Id, the effect of Ai on the gap is roughly additive with A2, while 
7] (dashed lines in Fig. Ic) tends to narrow the antibonding band and shift it upward (as 
found earlier [Q) without modifying the region near the vHs's. (To avoid an overall band 
shift, I have assumed A + 21] = 4eV for the dashed lines.) 

All of the above refer to the LTO phase. The effect of the Aj's on the LTT phase is nearly 
identical, as illustrated in Fig. 3a. The electron-phonon parameter rj has a much larger effect 
in the LTT phase Fig. 3b. For the present paper, however, my main interest in in the 
role of spin-orbit coupling on the LTO phase, so the LTT phase will not be discussed in as 
great detail. 
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E. Splitting of the vHs 



The nature of the sphtting of the vHs depends sensitively on the magnitude of the 0-0 
hopping parameter, too- This is illustrated for four different values of too in Figs. 4 (energy 
dispersion), 5 (density of states), 6 (Fermi surfaces) and 7 (electron filling). The peaks in 
the dos correspond to vHs's, at which the topology of the Fermi surface changes. There can 
be up to four of these peaks, at different energies, in the LTO phase. In this section, the 
origin of each peak will be discussed. Note also that for too < 0.1, the dos vanishes at a 
point between the vHs's. This is a signal for a diabolical point in the energy bands, which 
will be discussed in more detail below. 

It is easiest to begin with too—0, Figs. 4a, 5a, 6a. In this case, there is an approximate 
electron-hole symmetry, so doping only up to half filling need be discussed. For a nearly 
empty band, the first electrons go into pockets near the F and M points (all descriptions 
are based on the tetragonal zone - in the orthorhombic zone, these points are equivalent). 
With increased doping, these pockets fill up and expand toward the diagonal, Cx + Cy — 0, 
just as in the case without spin-orbit coupling. Near half filling, however, a first topology 
change occurs when small pockets appear near the X and Y points of the zone. This is 
shown in detail in Fig. 8a. At this point, all of the Fermi surfaces are closed and electron- 
like. Additional doping causes all surfaces to grow further, until the two surfaces meet to 
produce a rectangular box-like hole Fermi surface. This electron-hole crossover marks the 
point of the largest dos peak - the main vHs. With additional doping, there is a single 
hole-like Fermi surface (in this quadrant of the Brillouin zone), centered on the point M. 
As doping increases, the area of the Fermi surface shrinks, until at half filling it is reduced 
to a single point at M. This is the diabolical point: here, the spin-orbit coupling is unable 
to lift the degeneracy of the two bands, leading to a biconic structure, or 'diabolo', in the 
energy surface. At this point, the dos vanishes. For doping beyond half filling, the entire 
sequence is played out in reverse order. Fig. 7 illustrates the integrated electron density n 
directly. Since this is the integral of the dos, the steep rises correspond to the vHs peaks. 
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while the flat region is associated with the low dos regime near the diabolical point. 

For finite too, the doping dependence is similar, but more complicated, since there is no 
longer electron-hole symmetry: the diabolic point is shifted above the midway point between 
the two principal vHs's. Thus, in Figs. 6b-d, the topology changes are all similar for doping 
up to the lower main vHs, corresponding to the box-like or coffin-like Fermi surface (dotted 
lines) . This is presumably the Fermi surface of main interest for superconductivity, since it 
falls in the doping regime x ~ 0.02 — 0.25, depending on the value of too- Note also that 
increasing too decreases the aspect ratio of the box, leading to a somewhat squarer Fermi 
surface, with the corners moving away from the X and Y points (see also Fig. 8b). The 
long sides acquire some curvature, giving the Fermi surface a coffin-like appearance. 

New features arise in the more electron-doped regime. For too ==0.1, Figs. 4b, 5b, 6b, 
doping beyond the box vHs leads first to the diabolical point, and hence a zero in the dos, 
just as for too — 0. Beyond this point, however, a new topology change occurs, as electron- 
hke Fermi surface pockets open up near the X and Y points (dashed lines). Finally the 
main vHs on the electron-doping side arises when these corner pockets intersect the growing 
electron pocket centered on M (solid lines). 

Increasing too further, a point is reached when the electron pockets at X and Y appear 
before the diabolical point is reached (Fig. 6c, too = 0.25). Beyond this point, the Fermi 
surface becomes semi-metallic, with electron pockets near X and Y, and a hole pocket near 
M. The dos no longer vanishes, and half fiUing corresponds to the doping at which the hole 
and electron pockets have equal areas. Doping beyond half filling, the hole pocket shrinks 
to a point (the remnant of the diabolical point) (dot-dashed lines), then an electron pocket 
grows from M, until the electron pockets merge at the final vHs. The case for too — 0.65 is 
similar. Fig. 8b shows how the box- like Fermi surface evolves with increasing too- 

A similar series of topological changes arises in the LTT phase. Fig. 9. 
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III. DIABOLICAL POINTS OF ENERGY BANDS 



Whereas quantum mechanical energy levels tend to repel one another, to avoid degenera- 
cies, level crossings can appear for sufficiently complex Hamiltonians. The analysis of such 
level crossings in molecules is an important aspect of understanding reaction rates of chemi- 
cal transitions. For a diatomic atom the degeneracy of two electronic levels is extremely rare, 
but Teller showed that level crossings should regularly occur for larger, polyatomic molecules 
Tsj . Such crossings are closely related to the JT effect, and Herzberg and Longuet-Higgins 



demonstrated that the degeneracy leads to anomalous statistics: to make the electronic 
wave function single valued, the orbital angular momentum must be quantized in half- 
integral multiples of h. An alternative approach is to introduce a pseudo magnetic field, 
so that single-valuedness of the electron is guaranteed by associating a quantized vortex of 



pseudoflux at the molecular site [15 



The degeneracy of the atomic levels leads to a linear dispersion of energy away from the 
point of degeneracy, so that the energy surface has a biconical form - Berry has called such 
points diabolical points |]TB|]. A signature of such points is that the electronic phase changes 
by vr when the electronic orbit traverses a path in energy around the point [|r^. This phase 
change was an early example of what is now known as Berry's phase |T^. Diabolical points 



have also been found in generalized two-dimensional Hubbard models, in the presence of 



orbital antiferromagnetism or spin nematic order |T^, and similar phenomena are found in 
a number of elemental semiconductors, notably tin, which is a zero-gap semiconductor, and 
(approximately) graphite. A key parameter is the critical value of too, which the half 
filled band undergoes a semimetal (SM) - zero-band-gap-semiconductor (ZGS) transition. 
This will be calculated ClS db function of Aj below. 

The Fermi surfaces near a diabolical point have many of the right properties to explain 
a number of anomalous properties found in the cuprates, as summarized in the Discussion 
Section. For example, near the diabolical point, the LTO phase is stabilized by the umklapp 
splitting of the vHs; doping, if accomplished by a rigid band shift, leads to hole-pocket Fermi 
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surfaces, which provide the correct Hall density; and finally, sufficient doping repositions the 
Fermi level at the new vHs of the box/coffin-like Fermi surface, close to the doping of 
optimum Tc. 

The form of two-fold degeneracy illustrated in Figs, la, 2a, is not unique to the LTO 
phase of LSCO, but is a common feature of a number of space groups containing glide 
planes. In all of these cases, spin-orbit coupling lifts most of this degeneracy, leaving a 
residual two-fold degeneracy on a few isolated points or fines. In the case of LSCO, when 
the LTO energy bands are treated as two-dimensional, the residual degeneracy is a clear 
example of a diabolical point, as illustrated in Figs. 4a,b. The situation can be somewhat 
more complex in three dimensions. Thus, for LSCO, when c-axis dispersion is included, 
there is an isolated degeneracy (diabolical point) at one point along the line perpendicular 
to Y*, while the entire line passing through X* retains a two-fold degeneracy. Since this 
dispersion is small, its inclusion should not greatly change the results of the present paper. 

The diabolical point is a novel form of vHs, associated with overlapping bands. This can 
most easily be understood with reference to the dos. Fig. 5. In Fig. 5a,b the dos goes to zero 
at the diabofical point. This is the characteristic form of the dos at the top or bottom of a 
band, so the diabolical point can be understood as the point where one band terminates and 
another immediately begins. Note that the electron density n is defined as the number of 
electrons per Cu; since there are two Cu per unit cell in the LTO phase, the diabolical point 
appears when there are exactly two electrons per cell - an exactly filled band. However, 
when too increases, the topological change occurs while the two bands still overlap in other 
parts of the Brillouin zone. Figs. 5c,d. In this case, the diabolical point is shifted away 
from half filling. The local contribution to the dos from the region of the diabolical point 
vanishes, but there is a finite total dos due to the contribution from other parts of the Fermi 
surface. 



IV. SLAVE BOSON CALCULATIONS: CHOICE OF PARAMETERS 
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A. Spin Corrected Slave Boson Theory 



Slave boson calculations offer a systematic procedure for incorporating correlations into 
a tight-binding model of the Cu02 planes. The model has a number of desirable features, 
including renormalization of the bandwidth near half filling, and a metal-insulator (Mott) 
transition at exactly half filling above a critical value of the ratio t^jtcuo- A key limitation is 
that the theory is formally a large-N theory, where N is the degeneracy of the electron states 
on a Cu site, and the calculation is usually carried out only to lowest order (in 1/A^). At 
this order, the slave boson calculation does not include magnetic effects. Moreover, ~ 2 
for the cuprates, so it is doubtful whether the lowest-order predictions of the theory are 
quantitatively correct. Nevertheless, detailed calculations are important - for example, to 
see how close the cuprates are to the metal-insulator transition, for a 'realistic' choice of 
parameters. In the present context, it is also important to determine how close the state at 
half filling is to the diabolical point. 

In the large-N calculations, some parameters are renormalized by a factor A^, which 
then does not explicitly appear in the final Hamiltonian. Nevertheless, in choosing the 
'correct' parameters to use in the calculations, this factor should be included. In order to 
minimize confusion, I will write the 'large-N' parameter as primed, and the bare parameter 
as unprimed. In particular, in comparing the present results to e.g., Ref. one should 
take t'^^Q = ^J^tcuOi = 2V^, where V is the nearest-neighbor Coulomb repulsion. 

In conventional slave boson calculations, the factor A^ includes both spin and orbital 



degeneracy (for the cuprates, N = 2 due to spin degeneracy [^). However, I believe that 
it is not appropriate to include the spin degeneracy in the term t^^Q. The factor A^ arises 
as follows: if the electron can hop from the Cu to A^ different orbitals on a given O, and 
all of these orbitals have the same energy and hopping t, then the antibonding band of the 
A^ + 1 X A^ + 1 matrix has the same form as that of a 2 x 2 matrix, with a single effective 
O orbital, if the hopping energy is renormalized to V^^- Such an enhancement would arise 
in the cuprates only if spin flip hopping were as probable as spin-non-fiip hopping. But the 
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spin-flip hopping is considerably less probable, and in fact is explicitly calculated here as 
spin-orbit coupling. 

If spin-flip and non-spin-flip scattering are equally probable, then the electron loses spin 
memory on each hop, greatly modifying the magnetic properties of the model. Hence, I 
propose a 'spin-corrected' version of slave boson theory, in which spin-flip and non-spin- 
flip scattering are explicitly distinguished. In this case, the bare hopping parameter t'fj^Q 
becomes 

^ _ 

^cuo - y y^CwO- 

Since all bare parameters are normalized to tcuOi this modification does not directly affect 
any of the predictions of the model. However, it is crucial when one attempts to determine 
whether the Mott transition occurs in the cuprates, by estimating 'realistic' values for the 
parameters. For instance, the Mott transition occurs only if A is larger than a critical value, 
~ ?>.'ibt'Q^o- I^'o^ reahstic values of A, a transition will occur only if the factor \fN is not 
included. 

For convenience in the subsequent discussion, I will carefully distinguish the two cases, 
labelling the spin-corrected parameters with a superscript t (i.e., when it is assumed that 
^'cuO — ^Cuo), and without the superscript for the conventional slave boson theory (assuming 
t'^^Q = ^/2tcuo)- Thus, when Aq = Q\/2eV = S.SeV, then A^ = QeV . 

In contrast, the nearest-neighbor Coulomb repulsion V is the same for both spins, so the 
factor N = 2 should remain for V. 

In searching through the literature, I was unable to find any explicit discussion of the 
proper choice of degeneracy factors. However, it is clear that different groups have made 
different choices, and that this choice is a cause for significant variations in the resulting 
slave boson calculations. To be specific, let b be the slave boson amplitude, and ro the 
renormalization of the hopping parameter. Then b is defined through the equation 

N 

Y.dUi + b'^b^qoN, 
1=1 

where N is the electron degeneracy and go is the filling factor (for the present case, N — 2, 
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go = 1/2), while tq is defined by tn = r^tcuo- Given these definitions, then previous slave 



boson calculations have used a variety of values for the ratio r^b = ro/6, with Refs. [p^-|2^ 



assuming r/b = 1, Refs. [0,0 assuming r/b = 1/v^, and Refs. (apparently) assuming 

r/b = 1/2 (since Ref. uses the same symbol for tcuo and tj^^o, it is sometimes difficult 
to ascertain which is meant). In general, most of the calculations assuming 7]^ = 1 find that 
the metal-insulator transition does not occur at half filling, for reasonable parameter values, 
while the transition is found if a smaller 77^ is assumed. The present calculations confirm 
both these results. 

B. Bare Parameter Values: 'Large' Parameters 

While the technique of carrying out the slave boson calculations to lowest order in 1/A^ 
is reasonably well understood, and considerable progress has been made in calculating the 
next-order corrections, relatively little attention has been paid to the proper choice of the 
starting, bare parameters. This choice is crucial, particularly in the immediate vicinity of 
half filling, where it controls whether or not a metal-insulator transition will occur. Under 
the present circumstances, the introduction of additional parameters - to describe spin-orbit 
and electron-phonon coupling - renders a detailed discussion of the choice of parameters even 
more urgent. Two separate problems arise: how to choose the bare parameters, and how 
the parameters vary with doping. This section will attempt to provide an introduction to 
the problems involved, but cannot pretend to provide a definitive solution. 

The usual technique which has been adopted is to extract the bare parameter values 
from an LDA or cluster calculation of the band structure. However, simply parametrizing 
the the LDA-derived Fermi surfaces is not appropriate, since these calculations already 
include correlation effects in an average way. Fortunately, for LSCO a number of groups 
have analyzed the band structure calculations, to provide estimates for the bare parameters. 



Table I hsts the values produced by several groups, using LDA or cluster calculations p7|-p9 
Earlier calculations are summarized in Ref. IBIl . 
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Even knowing these values is not sufficient: it must be regognized that the three-band 
model is a simplification of the physics, and the parameters which enter it are effective 
parameters, which describe the combined effect of several microscopic parameters. Two 
examples are of particular importance. 

The nearest neighbor Coulomb repulsion V acts to renormalize the bare Cu-0 separation 



3l[, Ao, according to 

Ao^ Ao + 2V^(l-a;-4rg), (7) 

where tq is the mean field amplitude of the slave boson. While the term 2V{1 — x) is 
just a constant, and can be absorbed into an x-dependent Aq, the last term depends self- 
consistently on the renormalization, and hence gives rise to a possible instability against 



phase separation |^0|- However, the critical value of V required for this phase separation 



1.76t^„Q/2 ~ 1.6eV |^0| in the standard slave boson theory, or ~ m 



the spin-corrected version) appears to be somewhat larger than the value estimated from 
LDA-type analyses, and moreover near half filling the value of Tq tends to be very small. 
Hence, in the present calculations Eq. 7 will be approximated by enhancing the effective 
value of Aq. 

The second example involves the proper choice of the effective 0-0 hopping parameter. 
The LDA-derived value of too is too large to give a good description of the curvature of 
the Fermi surfaces, suggesting that higher-lying bands play an important role in modifying 
the Fermi surface shape Aligia has shown that this parameter is the sum of two 
effects of opposite sign: the direct 0-0 hopping, of value estimated in Table I, and the Cu 
dx^-y^ — interference term, which acts to reduce the curvature. 

The problem of calculating the effective value for too is similar to the problem of cal- 
culating t' m. a. t — t' — J model, since the role of both parameters is to distort the Fermi 



surface away from square at half ffiling. In this case, Jefferson, et al. have shown that 
even within a three-band model correlation effects can act to reduce t' toward zero, and can 
even change its sign. 
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In Appendix II, I derive an effective value of too from a five-band model of the Cu02 



planes, including the Cu (1^2 and apical O orbitals [Q. A substantial reduction is found, 
and too could even change sign. Due to uncertainties in many of the parameter values, 
it is difficult to use this calculation to fix the correct effective value for too- Hence, I 
propose an approximate criterion for an acceptible value of too'- that the vHs fall at the 
same doping as in the LDA calculations, for the self-consistently renormalized slave boson 
values. The rationale for this is that correlation effects will be smallest away from half filling, 
so that the approximate manner in which LDA calculations incorporate correlation effects 
are likely to be most nearly correct near the vHs. The remarkable degree to which these 
calculations reproduce the experimentally measured Fermi surfaces confirm the plausibility 
of this choice. This estimate of the effective value of too is in general agreement with the 
reduction estimated in Appendix II. 

Figure 10 shows how x^hs-, the doping at which the Fermi level coincides with the vHs, 
varies with too^ assuming bare values of tcuo = l-3el^ and A = 4 or 6eV. Using the bare 
parameters, the vHs falls in an acceptable regime (I estimate x^hs = 0.16 in LSCO and 
0.25 in YBCO). However, correlation effects greatly reduce the value of tcuO, and hence 
enhance the ratio too/tcuo, which controls the Fermi surface curvature. Thus, to recover 
reasonable values for x^hs, it is necessary to assume smaller values for too- For definiteness, 
I have assumed bare values of Aq = 4eV or 6eV, in which case ~ 0.14 — 0.2eV for 
LSCO, and 0.3-0.4eV for YBCO. 

It is tempting to speculate that the lower in LSCO is related to the finding that the 
effective too is only half as large in LSCO as in YBCO. A smaller too can reduce Tc directly, 
because the excess hole population is smaller at the vHs, or indirectly, since a smaller too 
means a more nearly square Fermi surface at half filling, thereby enhancing the possibility of 
a structural instability which competes with superconductivity. If too is lowered by repulsion 
from the d^^ level, then this parameter is particularly sensitive to the environment off of the 
Cu02 planes - e.g., other layers, and the apical oxygen. Indeed, there is a close sorrelation 
between A^, the d^^-yi — d^^ splitting, and the distance dA between the Cu and the apical 



oxygen (Appendix III). This correlation can be interpreted in terms of a static JT-like effect 
(N.B.: not to be confused with the vHs-JT effect) which produces the distortion of the six 
O's surrounding the Cu from a perfect octahedron [ p6| , p7| . Hence, anything which reduces 
the Cu-apical O separation also reduces the d — d splitting, in turn reducing the effective 
value of too- The correlation between small Tc and small optimum hole doping has been 
known for a long time, and it is also known that the apical oxygen plays some role, since a 
larger apical O potential correlates with higher Tc IS^ . 



C. Bare Parameter Values: Spin-Orbit and Phonon Coupling 

An important role in the theory is played by the magnitude of the spin-orbit coupling. 
Specifically, if A = Ai + A2 is large enough, the state at half filling will be a zero-gap 
semiconductor (ZGS) rather than a semimetal. The location of the zero-gap semiconductor 
- semimetal (ZS) phase boundary depends sensitively on the relative magnitudes of A and 
too- For A = 50meV, the ZS crossover occurs at too — 0.13el^; for A = 6meV, the crossover 
occurs for too about an order of magnitude smaller. 

The parameter A has been estimated by Bonesteel, et al. as 

A2 ~ {^)tcuoO, (8) 
9 

where 6 is the octahedral tilt angle, and Ag is the shift in the electronic (/-factor g due to 
spin-orbit coupling. With the estimates Ag/g ~ 0.1 tcuo = l-SeV (I use the bare 
value of tcuO to estimate the bare value of A), and experimental values for 9, Eq. 8 gives 
Xlto — GmeV, Xltt ~ 24meV . 

There are a number of factors which could substantially enhance the value of Xlto- First, 
if there is a dynamic JT effect, the experimentally observed value of 9 will be smaller than the 
microscopic local values, which would be more representative of local LTT-like distortions. 
Secondly, dynamical JT effects can greatly reduce the observed spin-orbit corrections to 
5^- values (the Ham effect [^). Thirdly, A is playing a role similar to umklapp scattering in 
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splitting the vHs degeneracy, and in the analogous one-dimensional problem, it is well known 
that renormalization effects greatly enhance an initially small umklapp term (for a discussion 
of the connection between the vHs problem and one-dimensional umklapp scattering, see 
pip. Finally, as will be seen in the following subsection, what is important in the present 
problem is how the parameters renormalize in the presence of strong correlation effects. I 
am assuming that A is renormalized in such a way that \/tcuO remains constant. This 
assumption is by no means self-evident - Bonesteel, et al. f^, suggest that A is unaffected 
by strong coupling. From the perspective of the present analysis, any increase in the ratio 
A/tc«o is equivalent to starting with a larger effective bare value of A. For all these reasons, 
it is plausible to assume that Xlto is large enough that the state at half filling is close to 
the ZGS phase. Clearly, a detailed microscopic estimation of the parameters, particularly A 
and too is a desideratum. 

The electron-phonon coupling parameter can be estimated as 

^ ^ 2^^^^ ^ 2.16^2e\/, (9) 

where ~ 1.6. This is close to the recent estimate of Song and Annett [|12|, rj ~ 2.W'^eV. 
Using Olto — 0.5 p|, tilto — ^AmeV. In the LTT phase, 9ltt is about four times larger, 
yielding t]ltt — 86meV. The latter value would also be approximately correct for a dynamic 
JT effect in the LTO phase. 



D. Doping Dependence: Scaling with tcuO 

In the absence of spin-orbit coupling and electron-phonon interaction, the slave boson 
theory provides a system of self-consistent equations for determining how A and tcuO vary 
with doping, due to correlation effects (Eqs. 11, below). An important question is how too, 
and, in an extended theory, A and rj vary with doping. This question can be conveniently 
separated into two parts: since strong correlation effects can renormalize tcuO to zero, while 
A remains finite, one must first ascertain how the parameters scale with tcuO, and then 
check if there are residual finite renormalizations, as with A. 
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The scaling can be determined from the theoretical expressions for the parameters, 

A ~ (10) 

and Eq. 6b for rj. Here Aq is the spin-orbit coupling constant |]^, is the hopping parameter 
between the planar O pa orbitals and the Cu d^^ orbital, and is the Cu (1^2 — d^^^yi 
splitting energy. Now correlation effects inhibit the placing of two holes on the same Cu, 
and this causes the renormalization of tcuo to zero at half filling. Thus, t' should also scale 
to zero, since it involves the same Cu orbital. The case with tz is slightly different, since 
it involves a different Cu orbital, but the Hubbard U for this orbital is also large, and the 
simplest assumption is that also scales to zero at half filling. 

Hence, the following assumption will be made for the scaling of the parameters. First, 
A is renormalized proportionally to tcuOi while r] scales as t^uO (0 Finally, too should not 
rescale with tcuO- 

Because of the strong rescaling of ?7, it might be thought that this electron-phonon 
coupling could be neglected near half filling, in comparison with the electron-phonon coupling 
term, A. However, as the scaling analysis of Section VB shows, the gaps produced by both 
terms scale in the same fashion near half-filling. Because of the different ways they enter the 
Hamiltonian, Eqs. 2 and 6, the effective electron-phonon interaction strengths are g ^ t] and 

~ Xtcuol ^1 both of which scale as g ^ near half filling. The electron-phonon coupling 
parameter has a similar scaling, Agp ~ g^N{Ep) ~ Tq, since N{Ep), the dos at the Fermi 
level, varies as N{Ef) ~ A similar scaling has been found in a number of earlier 



studies [42|. I believe that this suppression of electron-phonon coupling by correlations is an 



artifact of the mean- field theory. Anderson [43] has found that correlations actually enhance 



nesting effects, and Eliashberg |Q| finds that electron-phonon coupling can remain strong 
in the presence of strong correlation effects. This issue will be discussed further in a future 
article. 

A separate issue is the doping dependence of the LTO transition temperature, which 
drops to zero around x ~ 0.2. Since A ~ ^ and t] ~ 6"^, these parameters will suffer a further 
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reduction with doping. Ultimately, of course, the LTO transition temperature should be 
calculated self consistently in the presence of correlation effects. 



E. Residual Doping Dependence 

In addition to the strong rescaling noted above, there can also be finite renormalizations 
with doping, similar to the renormalization of A. As will be shown below, the renormalized 
A is typically a factor of three smaller than its bare value. This renormalization of A will 
have an impact on several of the other parameters. For instance, rj oc A^^ will be increased 
by the same factor of three. The effect on other parameters depends on how A^^ varies with 
doping. Since A^ << A, it is necessary for the d^^ level to shift in tandem with the dx2_y2 
level, to keep the states at the Fermi level of dx2_y2 symmetry. However, in this shift, it 
is not clear whether the two levels maintain a constant separation, or whether A^ is itself 
reduced proportional to A. The latter situation arises naturally when a single slave boson is 
introduced for all d-levels, and has been assumed in a recent study of the role of the apical 



oxygens |^ In this case, A oc A^ would be further increased by the same factor. 



A similar effect would arise for the effective value of too- While the microscopic too is 



unlikely to be strongly affected by Cu Coulomb correlations (but see [@]), the counterterm, 
due to repulsion from the dz2 levels, will be enhanced by reducing A^^. The form of the 
reduction is suggested in Table II (where an exlpicit form for X is given in Appendix II). 

In summary, there is much left to be understood about the proper choice of parameter 
values, and how they change in the presence of strong correlations. The fact that reduc- 
ing A increases some other parameters could even lead to an instability, just as occurs with 
sufficiently large values of V , Eq. 7. Lacking more precise knowledge, the following prescrip- 
tion will be followed here. Only renormalization of A and tcuo will be explicitly accounted 
for in the slave boson calculations. It will, however, be assumed that A scales with tcuO, 
and 7] with t^'no- This will be shown to lead to an interesting scaling regime when tcuO 
is sufficiently small. The other, minor renormalizations discussed in the present subsection 
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will be approximately accounted for by adjusting the values of the bare parameters. Since 
the renormalization is strongest near half filling, the bare parameters will be taken as those 
appropriate to half filling, and effects due to their variation with doping will be neglected. 
A summary of the present parameter estimates is given in Table II. 



V. SLAVE BOSON CALCULATIONS IN THE THREE BAND MODEL 

A. Slave Boson Technique 

The strong on-site Coulomb repulsion acts to greatly reduce the probability of double 
occupancy of an atomic orbital. The effect is most important for Cu, since U is larger for 
copper and since the Cu band is very close to half filled, so the probability of finding an 
unoccupied Cu site is small. Hence, for simplicity U will be neglected on the O sites, and 
set to infinity on Cu - this has the effect of forbidding any double occupancy of the Cu's. 
In a slave boson calculation, this constraint is enforced by introducing a boson to occupy 
the second orbital of each singly occupied site. By satisfying the constraint in mean field, 
it is possible to reduce the calculation to a self-consistent renormalization of the bands. 
Basically, the Cu-0 hopping parameter is reduced by the renormalization, while A tends to 
increase, as the antibonding band becomes more purely Cu-like. 

The equations of self-consistency can be written 

rl - 111 -E^lME,)], (11a) 

k 

Ao - A = -1^ E ulME,)iE, - A), (116) 

^^0 k 

where the renormalized parameters are tn — TqIq^q and A, fh{Ek) is the Fermi function, 
and Uk is the d-wave amplitude of the wave function. and Uk are calculated using the 
renormalized parameters, and it is assumed that too is not renormalized. Only the case 
T = will be considered in the present paper. As written, Eqs. 11 are valid for a hole 
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picture, so fh{Ek) = 1 for Ek > Ep, and = otherwise. There is an equivalent electron 
picture, discussed in Appendix IV. The sums are normalized per Cu atom, so for example 

E(l) = 6, (12a) 

k 

that is, that there are three bands (spin 1/2), and 

J2MEk) = l + x. (126) 

k 

Equations 11 are written in such a way as to remain valid in the presence of spin- 
orbit coupling and electron-phonon coupling. These equations are solved self-consistently 
by guessing initial values for the renormalized values and numerically performing the double 
integrals to allow calculation of the corresponding bare values, then readjusting the initial 
guesses until the known values of the bare parameters are recovered. [This improves upon 
my earlier work [HBl, which approximated the problem by a single integral.] If is small 



enough, however, the equations simplify to a scaling form, for which only a single parameter, 
A, need be varied. This scaling regime is discussed in the following Subsection. 

B. Scaling Regime 

When t/j ~ 0, and A >> 4too, the antibonding band is very narrow, 6 = E — A << E. 
In this case, the dispersion relations, such as Eq. 4, can be solved to lowest order in 6, 
setting E A in all the remaining terms. The resulting solution posesses a simple scaling 
form. 

It is convenient to first treat the special case, X = rj = - the untilted case. In this 
case, the normalized dispersion •y = 6/t (with t = 4t|,/A) depends only on the parameter 
y = Atoo/^- 

(13) 

The resulting bandwidth is 2t/(l — y), and the band is nearly pure Cu-like, with 

«^ = 1 - |, (14a) 
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Letting 7 = J2k lfh{Ek) be the average of 7 over the occupied states, with a similar definition 
of e, the self- consist ancy equations become 

Ao = A + %^, (15a) 

where x is the hole doping. Note that once A is known, Eq. 15b is an explicit equation for 
in, and that in the scaling regime t^j; ^ at half filling. 

Remarkably, using the scaling of the small parameters discussed above, the same scaling 
regime holds for the full dispersion relation, Eq. 4. Thus, Eq. Al in Appendix I is a 
quadratic equation for 5, and it can easily be seen that when A scales oc t/j, and 77 oc then 
5 oc t, as found above. Now, however, the energy dispersion is a function of three variables y, 
X/tji, and i]/t'j^. The bare parameter values are sufficiently close to this scaling regime that 
most of the results of Figs. 1-9 are only qualitatively changed if the renormalized parameters 
are used in place of the bare ones (except when too = 0.65e\^, for which Atpo << A is not 
always satisfied). 

C. Slave Boson Results: (a) Strong Correlation {t —>■ 0) Limit 

The above scaling results can be used to analyze the strong correlation limit - specifically, 
to explore the circumstances in which tji is renormalized to zero. In particular, Eq. 15b 
shows that tj? = can only occur at half filling, x = 0. This also follows directly from 
Eq. 11a when tq = and Uk = 1, and is consistent with the result for the one-dimensional 
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Equation 15a shows that there is a minimum value of Aq for which t^ = is possible. 
Defining A^ = 4t^„o7o, where 70 is 7 evaluated for a half filled band, then the minimum 
value of A is A^, and at this point Aq = 2Am- Numerically evaluating 70 as a function of 
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y thus provides a plot of Ao/t^no vs too/tcuo = V^ml^t'cuOi Fig- H- Fo^^ 1/ = 0, 7o can be 
evaluated analytically, yielding p5| 



Ao 1 2 , , 

4a/- + — = 3.353, (16a) 



t'cuO V 2 TT 



or Aq = 4.36ey when tcuo = l-3eV^. Note that for the conventional theory, Aq = V2 x 
4.36 = 6.17eV. If a larger value, tcuO = l-4eV^, is assumed, then Aq = G.GieV, in good 



agreement with the 6.5eV found by Sudb0 and Houghton [p^] . 

The above result is, however, only a lower limit to the critical Aq when too 7^ 0. The 
analysis assumes that the system is in the scaling regime, which is not necessarily the case 
when too is comparable to A/4. A more accurate phase boundary can be found from 
numerically solving the self-consistent Eqs. 11, and this is also illustrated in Fig. 11. Care 
must be taken in deriving these points, since the numerical calculations become unstable 
when tn ~ 0. The procedure followed here is to choose a value for the renormalized tn and 
adjust A until tcuO = l-3el^. A number of values of were chosen in the range 0. 075-0. 2eV, 
and it was found that a plot of Aq vs gave a straight line, which was extrapolated to 
tfi = to generate the points plotted in Fig. 11 as open circles. These closely approximate 
a straight line 

Ao = 3.353t;.„o + 2.94too. (166) 

This metal-insulator phase diagram is in semi-quantitative agreement with Monte Carlo 
calculations on the three band Hubbard model, at least at too = 0, where most of the 
calculations have been done. Assuming a large but finite value for the on-site Coulomb 
repulsion, Ud = QtcuO, a charge-transfer insulating gap appears to open up for A/tcuO 
between 1 and 4 or between 2 and 3 ||4^. The resulting gap values are also comparable. 



This good agreement holds for the spin- corrected theory, while the conventional slave boson 
theory is off by a factor ~ ^/2. 

The ZS crossover can also be determined analytically. The diabolical point is the solution 
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of Eq. Al (Appendix I) at the M point {cx = Cy = 0), or 

Ed5d = Yzr^2 ' 

where the subscript d refers to the diabohcal point, the various Aj's are defined in Appendix 
I, and y' = too /Ed- Equation 17a is true in general; in the scahng regime, each exphcit 
should be replaced by A. The ZGS state is present when there are no other solutions of Eq. 
Al along the diagonal Cx + Cy = 0. When such solutions exist, the diabolical point falls in 
a semimetallic regime. The crossover ooccurs when the second solution falls at the X ot Y 
point, in which case 

Ed5d = 2{Xl + 2t^) ± AV2t\. (176) 

The upper solution is relevant here; the lower one falls close to the box-like Fermi surfaces. 
Equating Eqs. 17a,b yields the equation of the crossover 

yX>^l + Xl + V2{X, + A2) + 1] + y[2 - A1A2] - 4y2(Ai + A2) = 0, (17c) 

where Aj = Xi/tcuo- The resulting crossover phase diagram is shown in Fig. 12. 

D. Slave Boson Results: (b) Metal-Insulator Transition 

The full Equations 11 have been solved for a variety of bare parameters Aq, tcuo and 
tooi as a function of hole doping x. The following results are found. (1) For x 7^ 0, the 
system is metallic, with > 0. (2) There is a critical value Ac (Eq. 16), such that t/j — 
as X — 0, if Aq > Ac. (3) For A > Ac, the Fermi level jumps discontinuously as x passes 
through zero, with the discontinuity being the charge-transfer gap. (4) The renormalization 
is such as to keep the Fermi level pinned near the vHs - the more so as A increases. This 
pinning was found in earlier slave boson p6| , p3| , |2^ and other calculations. 

Using the parameters estimated above (Table II), the present calculations predict a 
metal-insulator transition at half filling, but only if tCuO' = tcuO (the parameters with 
superscript t in Fig. 11). If, on the other hand, it is assumed that tCuO' = \^tcuO, then 
the half filled band remains metallic, A < Ac. 
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The metal-insulator transition is illustrated in Fig. 13, both in the presence of spin-orbit 
coupling, A2 = 0.05eV (symbols) and in its absence (lines), assuming Iqq = 0.25eV, as 
Aq increases from AeV (solid line or open squares) to 5eV (dashed line of filled circles) to 
6eV (dot-dashed line or open circles). The transition in the absence of spin-orbit coupling 
will be discussed first. Figure 13a shows that t/j renormalizes to zero exactly at half filling 
(Aq = 5eV is just below the transition, so has a small but finite value at half filling). The 
scaling theory, Eq. 15b, approximates the doping dependence of tpt near half filling (dotted 
line). Near this point, the Fermi energy changes rapidly (Fig. 13b), from a value near Aq 
(Cu-like carriers) when x < to a much smaller value (more 0-like) for x > 0. This change 
is discontinuous above the Mott transition; this is the charge-transfer insulator gap [^, and 



should be compared to the one- dimensional Hubbard model (see Fig. 3 of Ref. |^). For 
even larger x, A becomes negative, and the Fermi level falls directly in the 0-like band. [In 
this regime, the three-band model should not be trusted, because of the proximity of the d^^ 
band to the Fermi level. Indeed, in this doping regime an enhanced d^^ character is often 
observed experimentally.] 

The quasi-pinning of the Fermi level to the vHs is illustrated in Figs. 13c, d. In a rigid- 
band filling picture, the bandwidth and the position of the vHs would both be independent 
of doping, so the separation between the vHs and Fermi level would smoothly track the 
variation of the Fermi level with doping. This is not what happens in the presence of strong 
correlation effects. Figure 13c shows that the position of the vHs (defined by the hole 
doping x^Hs at which the vHs would coincide with the Fermi level) actually changes with 
doping, while Fig. 13d shows the energy separation between the Fermi level and the vHs, 
AE = Ep — EyHs- Note for example the case Aq = 6eV: over the doping range x = —0.1 
to a; = +0.3, the Fermi energy varies by 4eV, while |Aii^| < lOmeVl Thus, correlations 
strongly pin the vHs near the Fermi level, especially when a Mott transition occurs at half 
filling. 

Figure 14 illustrates the corresponding Fermi surfaces, both at the vHs (i.e., when AE = 
0) and at half filling. Note that these latter Fermi surfaces are far from square, and their 
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proximity to the vHs arises mainly from the very small bandwidth. 

The changes in tpi and Ep due to inclusion of spin-orbit coupling are relatively small, 
Fig. 13a,b. (It should be noted that the routine for solving Eqs. 11 does not converge well 
for large Aq near a; = 0, or for smaller Aq at large x.) There is a curious symmetry between 
the doping of the vHs in the absence of spin-orbit coupling x^hs and that of the diabolical 
point Xdia in its presence: Xdia — —x^hs, Fig. 13c. For the assumed value Iqq = 0.25eV, 
the system is a semimetal at half filling. Figure 15 shows how the bands and Fermi level 
evolve with doping for Aq = 6eV. At x = 0.3 (Fig. 15a), the Fermi level lies in the gap 
between the two vHs at the X-point, but there is still a cylindrical Fermi surface centered 
at M, dashed line in Fig. 15d. Reducing x to 0.15 (Fig. 15b) or 0.05 (Fig. 15c), the Fermi 
level has now crossed one of the vHs near X, leading to the presence of both electron and 
hole pockets. Fig. 15d. 

In Figure 16, the same sequence is repeated for tQQ = 0.1. As x ^ 0"*", too is still 
too large to have a zero-gap semiconductor, in agreement with the scaling results. Fig. 12. 
However, on the Cu side of the charge transfer insulator (x < 0), the much larger value of 
A leads to a smaller value of y, with a corresponding zero-gap semiconductor at x = 0~ - 
i.e., the diabolical point falls at the Fermi level at (slightly less than) half filling! Figure 16e 
shows the band structure for x = —0.05 - which is as close to half filling as the numerical 
routine can be successfully used. The diabolical point is clearly in a true gap in the dos. 
The Fermi surface for x = —0.05 is composed solely of an electron-like surface (dotted 
lines in Fig. 16f) - the lack of hole-like sections is further evidence that the net electron 
concentration at half filling can vanish only by the Fermi surface shrinking to a point (ZGS) 
and not by compensation of electron and hole sections (SM). 



E. Discussion: Mott Transition? 



Early slave-boson calculations ||2T|], which neglected direct 0-0 hopping {too = 0), identi- 



fied the Mott transition as the point at which tcuO is renormalized to zero. Later calculations 
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included direct oxygen-oxygen hopping, too HOSJH]) t)ut it was found that slave boson 



calculations lead to a metallic state at half filling when the bare value tpo = 0.65eV is 
assumed. As argued above, the value of too used in these calculations should be a smaller, 
effective value, incorporating the role of coupling to the ci^a states. Moreover, the value of 

in the relationship t^^o ~ VNtcuO should, as argued above be taken equal to 1, and 
not 2, as is usually done. When these two corrections are made, the present results give a 
metal-insulator transition at half filling, for reasonable parameter values. 

This result is consistent with the results of other, non-slave boson calculations of the 



three-band model. Thus, the studies that produced the parameters of Table I p7H30|] also 
explored the low energy sector of the three-band model, and all find an insulating state at 



half filling with a large charge transfer gap, comparable to experiment. See also Refs. [51 



which find similar results. For all of these parameter sets, the slave boson calculation with 
N = 2 would have incorrectly predicted a metallic state at half filling. 

Zhang, et al. recently suggested that the slave boson calculations are in error near 
half filling, and that the state at half filling is always insulating (at least if too is not too 
large). Furthermore, when the slave boson calculation does predict an insulating state at 
half filling, its detailed predictions are in good agreement with more accurate treatments. 
It is still possible that, for finite too, there is some minimum value Ac, below which there 
is no metal insulator transition, since the Bi and Tl based cuprate superconductors do not 
show evidence for an antiferromagnetic insulating phase. 



VI. CONCLUSIONS 

The present results have an important bearing on the question of the role of the vHs in 
driving structural transitions in the cuprate superconductors. A major stumbling block has 



been the assumption |jTO[ that the vHs degeneracy is not split in the LTO phase. Spin-orbit 
coupling does produce such a splitting, playing the role of the umklapp scattering postulated 
earlier [^. However, vHs's cannot be destroyed in a two-dimensional band, but only split 
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up. The ensuing Fermi surface reconstruction leads to a large variety of new vHs's. Perhaps 
the most interesting possibility is a diabolical point, at which the Fermi surface collapses to 
a single point. 

While these Fermi surface reconstructions are undoubtedly present, just how important 
they are depends on the proximity of the Fermi surfaces, as a function of doping, to the 
topological changes of the bands. The answer to this question, in turn, depends on the 
choice of the band parameters, and how they are modified by correlation effects. These 
are difficult questions, and the present work is intended as an introduction to the problems 
involved, as well as a first estimate for many of the parameters. 

Among the questions which depend on these parameter choices are the following. Will 
the state at half filling be metallic or insulating? Is the LTO transition driven by the vHs 
splitting? At what doping does the vHs coincide with the Fermi level. Figure 11 provides an 
answer to the first question: the critical value of Ao/t'c„o above which tcuo is renormalized 
to zero by correlation effects, at precisely half filling. The critical value appears to be in the 
range of parameter estimates for both LSCO and YBCO. The transition is not driven by 
Fermi surface nesting, in that the transition always appears at half filling, even though the 
nesting is better at other dopings. However, there is an element of nesting involved, in that 
the critical value of Aq increases rapidly as too increases from zero. 

While I have demonstrated that the vHs degeneracy is split by the LTO distortion, it 
remains to be seen whether the energy lowering associated with this splitting can drive 
the structural transition. This will be the case for one special doping, when the vHs lies 
exactly at the Fermi level in the undistorted structure. However, I suggested a more 
interesting possibility: that correlations can pin the Fermi level close to the vHs over an 
extended doping range, so that the LTO transition could be Fermi-surface-driven over this 
full range. While this pinning is clearly present. Fig. 13d, the Fermi surface at half filling is 
not close enough to the vHs (Fig. 14) for this simple model to work. However, the present 
calculations suggest two closely related possibilities. 

One possibility is that a diabolical point falls at half filling, so the half filled band is a 
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zero-gap semiconductor, and hence naturally an insulator at low temperatures. If doping 
then produces a simple rigid band filling, the resulting small hole pockets can explain the 
observed Hall densities, and yet if the gaps are small the Fermi surfaces will appear to 
be large. Sufficient hole doping would bring the holes to a new vHs, associated with the 
box/coffin-like Fermi surface. Fig. 8b. However, it seems that too is too large for the 
diabolical point to fall at x = 0"*". 

The second possibility, illustrated in Fig. 16, is that the diabolical point falls at x = 0~. 
In this case, the LTO structure would be stabilized both at x = 0~ and at x = x^hs, but not 
at intermediate doping. In such a situation, an electronic phase separation may arise, to keep 
part of the holes at the vHs density This possibility arises for realistic values of too, 
but requires a somewhat enhanced value for the spin-orbit coupling parameter A. However, 
in the analogous one- dimensional metal problem, it is known that umklapp scattering is 
strongly renormalized to drive the metal-insulator transition at half filling. 

A very attractive possibility is that the renormalization arises from antiferromagnetism. 
It is known [^,0] that antiferromagnetism also lifts the band degeneracy along the X — M 
face of the Brillouin zone. Attempts to self-consistently calculate antiferromagnetism - to 
see if the gap opening lowers the energy sufficiently to stabilize the antiferromagnetic state 
- generally find that the paramagnetic state is more stable. Part of the problem is that 
the gapping is imperfect: residual Fermi surface sections are left behind. If, because of spin 
orbit interaction, the system is already close to a complete gap (a diabolical point), then it is 
possible that the additional effect of antiferromagnetic correlations will open a complete gap, 
and thereby enhance the stability of the antiferromagnetic state. It should be recalled that, 
upon doping away from the insulating state, the first carriers produce small hole pockets 
near the M points, just as would be expected for the diabolical point In the problem of 
the quantum Heisenberg antiferromagnet, it is known that Berry phases lead to dangerously 



irrelevant couplings [55 
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VII. DISCUSSION 



Spin-orbit coupling introduces a major reconstruction of the Fermi surfaces of the 
cuprates, and hence opens up a wide variety of questions, as well as suggesting possible 
answers to a number of problems involving the cuprates. Here, I briefly indicate some of 
the major issues. It should be noted that incorporation of this scattering clearly demon- 
strates the important role of the vHs in the structural transitions, in addition to their role in 
superconductivity. In brief, the LTO (and LTT) phase(s) are the equivalent of the charge- 
density wave phases in the vHs problem, and the physics of the cuprates is dominated by 
the competition between superconductivity and density-wave formation, in close analogy to 
one-dimensional metals. 



A. Umklapp scattering and the Hall anomalies 

It has previously been suggested that umklapp scattering could provide an explanation 



for some of the anomalous features of the Hall effect in the cuprates [56|. Basically, by 
introducing gaps at the vHs's, this scattering would change the Fermi surfaces from large 
areas (proportional to l+x) to much smaller areas, ~ x/2. That spin-orbit coupling produces 
just such an effect is clear from, e.g.. Fig. 2, thus confirming the suggested explanation. 

Moreover, the temperature dependence of the Hall density can, in principle, be un- 
derstood as a measure of the temperature-dependent gap associated with the LTO phase 
transition. As T is lowered below the HTT-LTO phase transition, the gap opens gradually, 
converting the hole gas from the large-area Fermi surface to the hole pockets, and then, as 
T falls further, reducing the areas of the pockets. 
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B. Mott and Antiferromagnetic Transitions 



This gradual gap opening can explain many anomalous features of the normal phase 
above T^, and particularly in the lightly doped regime near x ~ 0. This would be particularly 
true if, as a result of the combined effect of spin-orbit coupling and large, Hubbard-like 
correlation effects, the Fermi surface at half filling falls at a diabolical point, with vanishing 
area (as in Fig. 4a,b). 

In particular, the mobility shows very little change at the LTO transition, but there is 
a metal-insulator transition (change of sign of the temperature coefficient of resistivity) at 
much lower temperatures. This can be understood as another manifestation of a slow gap 
opening (enhanced in this case by the fact that the gap vanishes at the diabolical points). 
Near the transition temperature, the gap affects only the immediate vicinity of the vHs's, 
but these holes make very little contribution to the mobility, because of the very strong 



inter- vHs scattering ||5^. Indeed, by eliminating the vHs as a source of scattering, this 
could actually enhance the mobility of holes on other parts of the Fermi surface. 

Gradually, however, the shrinking area of the residual Fermi surfaces with decreasing T 
would cause the mobility to decrease, ultimately leading to a localization transition at half 
filling, when the areas shrink to zero. 

Just as in the usual localization picture, the spins should freeze out prior to the final 
localization, leading to a local moment formation on the Cu's. This may explain why 
the Heisenberg model for the antiferromagnetic transition has been so successful in these 
materials. In the past, the success of this model has led to suggestions for 'two carrier 
dynamics' - in which the holes on the Cu were localized, while the doped holes on the 
oxygens formed a separate band. Here, it is seen that both features arise naturally from 
spin-orbit scattering in the LTO phase. 
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C. LTO-LTT Competition? 



Since spin-orbit coupling directly splits the vHs degeneracy, it is not necessary to pos- 
tulate the existence of a dynamic JT effect, in order for the LTO phase to be electronically 
driven. Nevertheless, frozen phonon calculations find that the LTT distortion is energeti- 
cally favored over the LTO or any other local configuration . Moreover, while spin-orbit 
scattering does not strongly affect the LTT- LTO splitting (Fig. 3a), the ordinary electron- 
phonon coupling provides an additional energy lowering for the LTT phase (Fig. 3b). Thus, 
the possibility of dynamic JT scattering remains, and will be reanalyzed once the dynamics 
of the transition to a uniform LTO phase is better understood. 

It must not be forgotten that even in molecular systems there is a strong competition 
between the JT effect and spin-orbit coupling, and it will require considerably more work 
to sort out their respective roles in the cuprates. 



D. Incommensurate Diffraction Peaks 

Neutron diffraction studies have found that in La2Cu04, the antiferromagnetic peaks are 
commensurate, at Q = {71/ a, it /a). As the material becomes Sr doped, the peaks split and 
become incommensurate, at Q ± 6{TT/a,0) or Q ± 5(0,7r/a), with the incommensurability 
6 increasing proportional to the doping x A number of attempts have been made to 
interpret this as due to a peak in the electronic susceptibility Xg associated with a nesting 
feature in the Fermi surface, calculated using the three band model, with too 7^ 0. 

There are certain difficulties associated with such interpretations, however. In the first 
place, in a weak coupling calculation the intensity is orders of magnitude larger for inter- 
vHs scattering (which in these models is commensurate at Q) than the nesting-associated 
scattering when the Fermi level is shifted off of the vHs (see Fig. 8 of [0]). In a marginal 



Fermi liquid model, this vHs intensity enhancement can be greatly reduced by quasiparticle 



lifetime effects p2|. In the doped material, the model predicts an incommensurate peak in 
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the susceptibility, in accord with experiment; however, the peak usually remains incommen- 
surate at half filling |^0[|. To correct this problem, it is necessary to assume ||3^ that the 
vHs's are very close to half filling (since the nesting feature that they produce is commen- 
surate), and this is inconsistent with LDA band structure calculations, if a semi-rigid band 
filling model is assumed. 

While the details of the doping dependence remain to be worked out, the present model 
offers an alternative possibility for explaining the incommensurate peaks: the vHs's are 
shifted off of the X and Y points, so that inter-vHs scattering is incommensurate, and the 
peak positions depend on the interplay of the correlation effects and tilting instability, which 
modify the shape of the Fermi surface. 

Here, I can only indicate that the observed Fermi surfaces are consistent with the direc- 
tion and magnitude of the distortion 6. From Fig. 17a, it can be seen that there are a large 
number (9) of inter-vHs scattering vectors near Q, of which two are illustrated. Figure 17b 
shows the relative distribution of the nine peaks about the central commensurate position. 
The four outermost points match the positions of the observed incommensurate neutron 
scattering peaks (Fig. 17c). As far as the magnitude of the effect goes, the observation that 
5 ~ 2x would require the vHs peak to fall a fraction x of the distance from, e.g., Y 
toward F. For x = 0.14, this is about twice as far as the largest Fermi surface of Fig. 8b. 

While this result is suggestive, a number of problems remain. A more detailed calculation 
is required to determine the relative intensities of the various peaks. Moreover, it is not clear 
how the doping dependence of the peak arises in the present model. 



E. Phase Separation 

An additional complication which has been neglected in this paper is the possibility of 
some form of hole phase separation away from half filling, due either to magnetic effects ||61|| , 



electron-phonon interaction [^, or long-range Coulomb interaction (specifically, the next- 



nearest-neighbor Coulomb repulsion, V) |2^. In this transition, the holes bunch up in such 
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a way that part of them remain pinned at the insulating phase at half filling. The entries in 
Table II suggest yet another form of instability, which may be called a 'molecular JT' effect, 
since it arises when A or A^^ —>■ 0, which leads to additional electronic degeneracies. There is 



considerable experimental evidence for such phase separation ||6^,^J6^. If a self-consistent 
theory for the LTO structural transition can be developed, it will be possible to calculate 
the criteria for many of these instabilities. 



The physics underlying the electron-phonon induced phase separation |^ can readily 
be understood. The structural instability is stabilized by lowering the kinetic energy of the 
occupied electronic states, and hence the energy lowering is optimized when the Fermi level 
falls precisely at the (undistorted) vHs - or equivalently, half way between the two vHs of 
the distorted structure. In this case, the parameter regime associated with Fig. 16 is most 
interesting. There are two dopings at which the Fermi level falls at an optimal position for 
structural instability - both at x ^ 0.15, the original vHs, and at x = 0~, the diabolical 
point. The new point, a; = 0^ will be associated with the greater structural instability 
(higher transition temperature), since all of the dos is shifted away from the Fermi level. 
In contrast, for the old vHs, only part of the Fermi surface is gapped, leaving behind hole 
pockets. For intermediate dopings, the dos (in the distorted phase) increases, leading to a 
reduced kinetic energy lowering. If this reduction is large enough, the system will prefer to 
phase separate, keeping some holes at x ~ 0~ and the rest at x ^ 0.15. 



F. Future Directions 

Clearly, much work remains to be done, including working through detailed calculations 
of the effects discussed in the earlier subsections of this chapter. It will also be important to 
study the temperature and doping dependence of the LTO phase transition, in the presence 
of correlation effects. Fitting the positions of the incommensurate neutron peaks may allow 
some of the parameter values to be pinned down. It would also be desirable to extend the 
model to YBa2Cu307_5, where such incommensurate modulations have not been observed. 
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Appendix I: Band Dispersion Relations 

While the eigenvalue equation, Eq. 4, must be solved numerically, it is convenient to 
rewrite it as a polynomial in E. First, neglecting electron-phonon coupling (r; = 0), 

4E^Ai + 4ExA2 + x^As + 2(A - E)Bi + (A - EyB2 = 0, (^1) 

where x = 4:too, Si = sin{kia), and 

Al = X\ — (|Aa;|^Ca; + |AyPCy)^ + 2t^(2|Aa; — Xy\^ — X\)CxCy 

+2f{xi + \xi\-si + \Xy\x) + 1'(4 - (c. + -cyn m 

A2 = XU\X,\'sl + IXyl'sl) + t\\X. + Xy\' + 2t')isl + si) 

-tMXx\'-\Xy\'){c^^-cl), {A3) 

As = slsl{Xl + 2f)\ (AA) 

Bi = 2E'^{X\ + 2f) + E^x{Xl + 2t^){l + c^Cy) 

-Eix'mX^lX + \Xy\'sl + t'[sl + si]) - {xyA)slsliXl + 2t% (A5) 

B2 = E^- E\xy2){l + c^cy) + x%sl/16, {A6) 

and XI = |A,|2 + |AJ2, A| = A^A; + A^A^. 

Including electron-phonon coupling, the following terms should be added to Eq. Al: 

rJ'C^-7^''B2, {A7) 

where 

Cl = 4.E^(f[Cx + Cy] - \Xx\% - \Xy\%) + AE^Xt^iCx + Cy) 
+ Ex'^{f[Cx + Cy][l - C^Cy] - slCylX^l"^ " ^CxlAj/H - (XV8)A|(C^ + Cy), {A8) 



38 



and Tj' = 2rj{cx + Cy) in the LTO phase, or rj' = 2'qcx in the LTT phase, with ?7ltt = '^Vlto- It 
should be noted that, despite their comphcated form, Eqs. A4,A7 are a quadratic equation 
in Cy{E, Cx), except in the LTO phase when r] ^ 0. 

These equations also have an important scaling form, in the narrow band limit. That is, 
if t is so small that E » E — A, then E may be replaced by the fixed value A everywhere 
in Eqs A4-A9, except when it occurs in the combination E — A. In this case, scaling 
t ^ t' = rot, \i ToXi, for all i, r] ^ rlrj, and {E — A) — > rl{E — A), the system of 
equations is invariant. That is, in the scaling limit, only the bandwidth is rescaled, but 
the shapes of the Fermi surfaces are unchanged, and the dos keeps the same shape, but is 
enhanced by a factor l/r^. 

Appendix II: Effective Value of too 

Ila: Equation for tQQ 



Feiner, et al. introduced a five band model of the Cu02 plane, to explore the role of 
the apical oxygen. Here, I use the corresponding one-electron band structure to demonstrate 
the reduction in value of the effective 0-0 hopping parameter, too- The band dispersion 
can be written in the form 

A{1 - b[sl + si] - cslsD = 0, (51a) 

with 

A = E'^{E - A)V, (Bib) 
Ah = E[{E - A)W + AEtl^oVl {Blc) 
Ac = 9.[too{E - A) + 2tl^oWtooV - W], (Bid) 
V={E-A,){E-6o)-al 
W = 4[{E - A,)al + {E- 5o)al - 2aia2as]. 
The definition of the parameters differs somewhat from Feiner, et al. [BH]. The zero of 



energy is taken at the planar O p-levels, and too is taken as a positive number. The new 
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parameters are = A — A^, the Cu 6/^2 level, Sq, the apical O level, tcuA = "\/3tti, the Cu 
- apical O hopping parameter, to a = ^2/2, the apical O - planar O hopping parameter, 
and tz = a^i the Cu d^^ - planar O hopping parameter. Molecular orbital theory would give 
tz = tcuo I "\/3- Calculated values for these parameters are listed in Table 1. 

From Eq. Bla, the vHs is determined as that energy for which the point (e.g.) s^, = 1, 
= is at the Fermi level, or 6 = 1. For this choice of b, the shape of the Fermi surface, 
ky{kx) is given by the solution of Eq. Bla. Hence, all models with the same value of c will 
have the same shape of Fermi surface at the vHs\ In particular, the vHs will fall at exactly 
the same doping. Thus, for a three-band model to agree with the five-band model, the 
parameter too must take different values in the two models. [More generally, the three-band 
too will have a different value from the LDA too ^ presumably, the five-band model will be 
closer to the LDA value.] 

In the three-band model, the band dispersion is given by Eq. Bla, with 

c = 2y + y"", 

and y = AtQQ/EyHs- Hence the vHs of the three-band model will fall at the same doping as 
that of the five-band model, if too chosen to satisfy 

Khs^ 

where y' = Atoo/E'^us the five-band model. Note that not only is too 7^ ^c/o' t)ut the 
vHs is located at a different energy. In the five-band model, the vHs is given by the largest 
zero of 6 = 1, or 

{E[E - A] - Atl^o)y - - A)W = 0, (S3) 

while in the three-band model, the vHs is given by Eq. B3 with W = V = 1. Defining 
rjvHs = EyHs/ E'j^Hsi^ 1)) then Eq. B2 can be put into a suggestive form when ai = 0: 

This is essentially the form suggested in Table II. 

40 



Note that the above reasoning imphcitly assumes that the only role of too in the theory 
is to move the vHs away from half filling. This is in fact my present belief. If it is found 
that too plays some other role (or that, e.g., the apical oxygen plays such a role), then this 
would be evidence that the simple three-band model cannot capture the essential physics of 
the process. Until such time as that becomes necessary, however, it is simpler to work with 
the three-band model, with correctly chosen effective parameters. 

lib: Numerical Estimates of too 

Given values for all the band parameters, tQQ can be calculated from Eq. B2. Once 
more, there are separate problems in estimating the bare parameters and the renormalized 
parameters. Here, I will only give an illustration of the orders of magnitude expected. The 
new parameters will be assumed to have the values given in Table I [Q, with A = AeV, 



Az = 3.36eV", tcuo = 1-SeV, and too = 0.65eV^. Since the effect of each of the three 
hopping terms is quite different, I first illustrate how too would change if only one of the 
a's were non-zero (with value given by Table I). If only ai is non-zero, it has a small effect 
on the vHs, too ~ 0.57eV^. Either 02 or have comparable effects, with either leading to 
fJl = OAOeV. For all three a's non-zero, too ~ 0.28eV. Hence, the expected reductions 
can be quite substantial. For somewhat larger a values, too could even change sign. Thus, 
for 0^2 = 1, as = 0.75 (the molecular orbital value), and ai = 0, too = —0.1. 

Correlation effects will modify this result, particularly near half filling. In particular, 
both «! and should scale to zero as ~ tq, but a2, which involves an 0-0 hopping, should 
not be much affected. Since A is itself reduced by a factor of ~ 3, and it enters into the 
denominator of Eq. B2 or B4 (E' ~ A), then too niay not be significantly renormalized. 

In principle, going from a five-band to a three-band model could also renormalize the 
total width of the band - i.e., lead to an effective value for tcuo as well. In practice, this 
does not seem to occur, so only too is renormalized. 

Appendix III: Distortion of CuOe Octahedra 
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Ohta, et al. p8[, by compiling data on a large number of cuprate superconductors, were 
able to find only one consistent correlation between and a molecular property: the larger 
the apical oxygen - planar oxygen energy splitting, AVa, the higher is Tc. In the text, 
I suggest a related correlation: a smaller value of the d — d splitting AVd should cause a 



smaller effective value of too PS], which in turn correlates with a lower value of T^. 

The correlation of AV^ with da has an origin reminiscent of a 'classical' JT effect, and 
indeed Cu is a JT ion, often found in a distorted octahedral environment . That is, if the 



six oxygens around each Cu formed a perfect octahedron, then the Cu d^^-yi and dz^ levels 
would be exactly degenerate. The structural distortion lifts the degeneracy, as in a static 
JT effect. However, in the ordinary JT effect it is the splitting of the electronic degeneracy 
which is the driving force for the distortion. Here the problem is more complex - there are 
large commensurability strains between the Cu planes and the La planes, and the octahedral 
distortion leads to a considerable strain relief (this is why all octahedra distort along the 
same direction). The same sort of distorted octahedra appear in LaNi04, and Ni is not a 
JT ion P7|. 



Nevertheless, the distortion does cause the d-d splitting, as can be seen from an analysis 
of the data base of Ohta, et al. In particular, AV^ is linearly proportional to dA/dp, the 
octahedral distortion {dp is the planar Cu-0 distance). Figure 18a. Most of the compounds 
analyzed have only a five-fold oxygen coordination about the Cu, so even when dA/dp = 1, 
the two d levels are not degenerate. That is why the solid line in Fig. 18a does not pass 
through AVd = when dA = dp. For the 6-fold coordinated materials (circled in Fig. 18a), 
it is possible to draw an alternative line (dashed line) which does pass through this point. 

Now, when dA 7^ dp, the apical and planar O's will have different Madelung energies, 
so there should be a correlation between AVa and dA/dp. Unfortunately, this correlation is 
considerably weaker (Fig. 18b) than that with AV^, presumably because the apical O's see 
the non-Cu02 planes, which vary considerably from one material to another. 

Appendix IV: Electron Picture of Correlation Effects 
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Equations 11, which are vahd in a hole picture, can be replaced by the following equations, 
in an electron picture: 

rl^liE^lfeiE,)-!], {Dla) 

k 

^0 - A = ^ E«I/e(^^)(A - E,). {Dlb) 

k 

The equivalence of the two pictures can be shown by subtracting the corresponding pairs of 
equations (e.g., Eq. 11a and Dla) and using the fact that /e(-Efc) = 1 — fh{Ek) (valid for all 
temperatures). This results in 

^ k 

A^W4Ek. {D2h) 
^ k 

In the above equations, it is essential to keep in mind that the summation is over all 
fc-states of all three hybridized hands, Eq. 12. Equations D2 have a simple interpretation, 
saying that the three bands contain exactly one Cu electron (and two O electrons) per unit 
cell (Eq. D2a), and that the average Cu energy is unchanged by hybridization (Eq. D2b). To 
prove Eq. D2b, the explicit form of U]~ and E]^ are required. For simplicity, only the simplest 
case, based on Eq. 1 with too (and U) set equal to zero. In this case, the nonbonding band 
is purely 0-like, and may be neglected. For the bonding and antibonding bands. 



E,^^±J{^Y + Ul^olk, iD3a) 



7fe = sin^{kxa/2) + sin^{kya/2), and 



Using Eqs. D3, 



k 



^ 2Ek-A 



= [A + |(2 + 0)] = 2A, 
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thus confirming Eq. D2b. Here the prime means summation over two bands only. 

In the above argument, it is essential to note that the equivalence works only if the 
summation is over all three bands (in the more general case too 7^ 0). Hence, at T ^ 
the hole picture is considerably more convenient, since the summation can in this case be 
restricted to the antibonding band only. 

Table I: Bare Parameters (eV) 

Parameter Reference 



mm m 



A 


3.6 


3.51 


3.5 


tcuO 


1.3 


1.47 


1.3 


too 


0.65 


0.61 


0.65 


V 


1.2 


0.52 


< 1 


A, 




0.64 




5o 




1.46 




«! 




0.47 








0.66 




as 




0.50 
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Table II: Suggested Effective Parameter Values 



Parameter 


Bare Value 
(eV) 


Scaling* 


Corrected Bare Value 
(eV) 


tCuO 


1.3 


~ VN'ro (Eq. 11a) 




A 


Ao + 2V 


(Eq. lib) 


5-6 


too 


tl,o ^ 0.65 


tho - X/A 


0.14 (LSCO) 
0.3 (YBCO) 


V 


0.005 (LTO) 
0.08 (LTT) 


~ N'rieyA 


0.016 (LTO) 
0.25 (LTT) 


A 


0.006 (LTO) 

0.021 (LTT) 


- ^fWr^e/Aa 


? 



*: N' — N/2 in the spin-corrected model. 
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Figure Captions 

Fig. 1 Energy dispersions of the LTO phase of the 2D Cu02 sheets, for various choices 
of band parameters. [Dot-dashed lines = dispersion along the line X —>■ M.] Parameters 
are t — l.SeV, A = AeV, too — 0.65eV in all cases, with A2, Ai, 77 = (a) 0, 0, 0; (b) 6meV, 
0, 0; (c) sohd hnes = 50meV, 0, 0; dashed lines = 50meV, 0, lOOmeV; (d) 50meV, 50meV, 
0. 

Fig. 2 Fermi surfaces of the LTO phase of the 2D CUO2 sheets, assuming the same 
parameters as Fig. la,b, respectively. The Fermi energies are (following the solid surfaces 
of Fig. 2a, from lower left to upper right) Ep — 4.5, 5.0, 5.28, and 5.5eV. 

Fig. 3 Comparing the Fermi surfaces of the LTT (sohd and dashed lines) and the LTO 
phases (dotted lines) [dashed lines = dispersion near Y; solid lines = near X]. (a) = same 
parameters as Fig. Ic (sohd lines); (b) = same parameters as Fig. Ic (dashed lines). 

Fig. 4 Energy dispersion in the LTO phase, for several values of too- Parameters are 
the same as the solid lines in Fig. Ic, except too — (a) 0, (b) 0.1, (c) 0.25, and (d) 0.65eV 
[N.B.: Fig. 4d is identical to Fig. Ic, sohd hnes]. 

Fig. 5 Density of states (plotted against electron density) for the same four too values 
as in Fig. 4. 

Fig. 6 Fermi surfaces associated with changes of topology (vHs or diabolical point), for 
the same parameters as Figs. 4,5. (a) Ep = 5.225 or 5.337eV (solid line - the hole-like and 
electron-hke Fermi surfaces have identical shape), E^ia — 5.2826^; (b) 5.226 (dotted line), 
5.338 (dashed), and 5.34eV (solid), Edia = 5.322eV; (c) 5.227 (dotted), 5.337 (dashed), 5.386 
(dot-dashed), and 5.393eV (sohd); (d) 5.230 (dotted), 5.337 (dashed), 5.581 (dot-dashed), 
and 5.583eV (sohd). The dot-dashed hnes are the diabohcal points. 

Fig. 7 Carrier density n (electrons per Cu) for the same parameters as in Figs. 4-6: 
from left to right, too = 0, 0.1, 0.25, and 0.65eV. 
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Fig. 8 Evolution of Fermi surfaces near primary hole vHs. (a) Parameters of Fig. 6d 
{too = 0.65eV), for several energy values, Ep = 5.21 (dotted line), 5.228 (dot-dash), 5.23 
(solid), and 5.25eV (dashed), (b) vHs for the four cases of Fig. 6: too = (solid line), 0.1 
(dashed), 0.25 (dot-dashed), and 0.65eV (dotted). 

Fig. 9 Fermi surfaces associated with changes of topology in the LTT phase, for the 
same parameters as Fig. 3b. Energies are Ep — 4.96 (dashed hues), 5.35 (dotted), 5.47 
(solid), and 5.58eV (dot-dashed). 

Fig. 10 Doping x^hs at which the Fermi level coincides with the vHs, vs too, assuming 
tcuO — l-3ey. Bare parameters = solid (A = AeV) and dashed (A = 6eV) hues; renor- 
malized parameters = open (Aq = 4:\/2eV) and filled (Aq = 6\/2eV) squares; alternatively 
normalized parameters (too) — open (Aq = AeV) and filled (Aq — 6eV) circles. Dashed 
hues = guides to the eye. Solid horizontal hues = expected values of x^hs for LSCO and 
YBCO. 

Fig. 11 Threshold of metal-insulator transition. Open circles = full slave boson cal- 
culation; hues = scahng regime. Parameters (A2,77o) = (0,0) (sohd line and open circles); 
(0.005eV,0) (dashed line); (0.05eV,0) (dot-dashed line); (0.05eV,0.1eV) (dotted line). 

Fig. 12 Semimetal (SM) zero-gap semiconductor (ZGS) crossover, in scaling regime, 
assuming Ai = (solid line) or Ai =A2 (dashed line). In the scaling regime, y = Stoo/Ao- 

Fig. 13 Slave boson calculation of metal-insulator transition, assuming Aq = 4 (solid 
line), 5 (dashed line) or 6eV (dot-dashed line): (a) tu; (b) Ep; (c) x-vHs (for the hues) or 
^dia (foi^ the symbols); (d) AE = Ep — E^hs (lines) or Ep — E^ia (symbols). The dotted 
line in Fig. 13a is the scaling regime approximation. The vHs is at the Fermi level when 
the corresponding hne crosses the dotted line in Fig. 13c. 

Fig. 14 Fermi surfaces corresponding to Aq = QeV in Fig. 13 (dot-dashed lines), for x 
— XyHs (dashed line), 0+ (solid hne), or 0~ (dot-dashed hne). 

Fig. 15 Renormalized band structure and Fermi surfaces corresponding to Aq = 6eV, 
too — 0.25e\^, A2 = 0.05e\^, for various dopings x — 0.3 (a), 0.15 (b), or 0.05 (c). Horizontal 
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line = Fermi level. Fig. 15d = corresponding Fermi surfaces: x = 0.3 (dashed line), 0.15 
(dot-dashed), 0.05 (solid). 

Fig. 16 Renormalized band structure and Fermi surfaces corresponding to Aq = 6eV, 
t^oo = 0.1e\/, A*2 = 0.05e\/, for various dopings x = 0.25 (a), 0.15 (b), 0.05 (c), 0.02 (d), 
r -0.05 (e). Solid horizontal line = Fermi level; dashed line in Fig. 16e = diabolical point. 
Fig. 16f = corresponding Fermi surfaces: x = 0.25 (dashed line), 0.15 (dot-dashed), 0.05 
(solid), -0.05 (dotted). 

Fig. 17 (a) Schematic of full Brillouin zone, corresponding to Fig. 8b, illustrating some 
possible incommensurate, inter-vHs scattering vectors, (b) Corresponding distribution of all 
possible scattering vectors around central commensurate peak Q = (vr/a, vr/a). (c) Observed 
incommensurate diffraction pattern in LSCO + indicates position of Q. 

Fig. 18 Energy level splittings, either Cu (1^2 — d^^-yi (a) or apical vs. planar O (b), 
plotted against normalized Cu-apical O distance, dA/dp. Letters refer to different super- 



conducting compounds, after Ohta, et al. [3S]. Circled letters are compounds with 6- fold 
(distorted octahedral) coordination of the Cu. 
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